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Abstract 

A new nonlinear electromagnetic wave mode in a magnetized dusty plasma is predicted. Its 
existence depends on the interaction of an intense circularly polarized electromagnetic wave with 
a dusty plasma, where quantum electrodynamical photon- photon scattering is taken into account. 
Specifically, we consider a dusty electron-positron-ion plasma, and show that the propagation of 
the new mode is admitted. It could be of significance for the physics of supernova remnants and 
in neutron star formation. 
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I. INTRODUCTION 
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Dusty plasmas play an important role in planetary as well as in astrophysical systems 
However, many prominent astrophysical environments in which dusty plasmas 
can be found are highly energetic and strongly magnetized, such as supernovae, neutron 
stars, and young stellar objects (YSO). There are important new wave modes is such dusty 
plasmas jfj. and the dust grain speeds may even become relativistic Jd[ |b|. Supernovae, in 
which extreme energy scales are reached, allow for many exotic phenomena. In these violent 
events, neutron stars, YSO and pulsars are formed. Quantum electrodynamics (QED) is 
an indispensable explanatory model for much of the observed phenomena around these 
highly condensed objects. Scattering of photons off photons is predicted by QED, and it 
can be a prominent component of pulsar physics, since pulsars offer the necessary energy 
scales for such scattering to occur. Related to the scattering of photons is the concept of 
photon splitting in strong magnetic fields [7J. It has been suggested that this effect could 
be important in explaining the radio silence of magnetars Moreover, the propagation 


of intense electromagnetic waves in plasmas gives rise to many interesting phenomena Jffl. 
The existence of new plasma modes due to QED effects has been noted recently [111 I12L Il3|. 
In the present Letter, we point out the existence of a different electromagnetic wave that 
may exist due to the interaction of charged dust and photons with the nonlinear quantum 


vacuum. 


II. BASIC EQUATIONS 


First, let us consider the classical propagation of a circularly polarized electromagnetic 
wave in a multi-component plasma embedded in an external magnetic field B 0 z, where B 0 
is the strength of the external magnetic field and z is the unit vector along the z axis. The 
dispersion relation for a circularly polarized wave with frequency u and wavevector kz is 
then given by (see e.g. Refs, jlj] and Q) 
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where the sum is over all particle species j, u ; CJ - = qjB 0 /m 0 j and uj P3 = {n^q 2 / eonn^) 1 / 2 is 
the gyrofrequency and plasma frequency, respectively, 'yj is the relativistic gamma factor of 
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species j, and the upper (lower) sign in the denominator corresponds to left (right) hand 
polarization. Here, noj denotes particle density in the laboratory frame, qj charge, m,Qj 
particle rest mass, and c the speed of light in vacuum. The plasma is assumed to consist 
of an admixture of electrons (e), positrons (p), ions (*), and charged dust particles (d). 
We note that the dispersion relation (JTJ) has been derived without applying the usual small 
amplitude linearization. For a given k, the polarization is determined by the sign of aIn 
the background state, the plasma is close to equilibrium, and we have approximate charge 
neutrality, J2jQj n 0j ~ 0- Thus, for strongly magnetized relativistic electrons, positrons and 
ions, a low-frequency (a; -C kc and cny eiP) j <C uJ ce , cp , c i) wave will satisfy the dispersion relation 
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which is appropriate for large amplitude electromagnetic waves in supernova. Here, 7 ^ = 

(1 + z/ 2 ) 1 / 2 , where 
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gives the relativistic gamma factor for dust with charge q d = Z d e, where e is the magnitude 
of the electron charge. Henceforth we will allow u < 0, thus incorporating the left hand 
polarization as negative frequencies, and removing the =|= in Eq. Q. 

Next, we investigate the effects of the nonlinear quantum vacuum on the dispersion 
relation @ of a circularly polarized electromagnetic wave in a dusty magnetoplasma. 

The weak held theory of photon-photon scattering can be formulated in terms of the 
effective Lagrangian density 

i? = ^ 0 + ^HE, (4) 


where «5f 0 = — tCq F a bF ab = |eo(E 2 — c 2 B 2 ) is the classical free held Lagrangian, and 
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is the Heisenberg-Euler correction 3,ll3], where F ab = \e abcd F cd , and \F ab F ab = —cE • B. 
Here, k = 2a 2 /r 3 /45mQ e c 5 ~ 1.63 x 10 _ 30 ms 2 /kg, a is the hne-structure constant, and h is 
the Planck constant. With F ab = d a A b — d b A a , A b being the four-potential, we obtain (see, 


e.g. Ref 
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where j a is the four current. 

For a magnetic field-aligned circularly polarized electromagnetic wave with the electric 
held E 0 = E 0 (x + iy) exp [ikz — iut) , the invariants satisfy F cd F cd = —2Eq (1 — k 2 c 2 /to 2 ) + 
2c 2 Bq and F cd F cd = 0. Thus, Eq. © can be written as 
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in the Lorentz gauge, and □ = d a d a . For right-hand circularly polarized electromagnetic 
waves propagating in a magnetized cold dusty pair-ion plasma, the matter effects can be 
absorbed in the wave operator on the left-hand side by the replacement □ —» —D(u;,k), 
where D is the plasma dispersion function, given by (see, e.g. Refs. 0,13, Q], and Eq. ©) 
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We stress that the application of Eq. © requires held strengths smaller than the Schwinger 
critical held E s = ml e c 3 /eh ~ 10 18 V/m [see, e.g. Ref. JljJ). Looking for low frequency 
modes, the dispersion relation obtained from Eqs. © and © then reads 
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Equation a is the main result of our calculations. Due to their large mass, the dust particles 
can under most circumstances be considered non-relativistic, i.e. ^ d ~ 1 ja, 6]. We note the 
following for the dispersion relation ©: In the classical case, according to Eq. m we must 
have u'yd ~ wcd < 0. This requirement is obviously not necessary for the QED case, where 
u'fd — io cd > 0 is allowed. Thus, new low frequency waves are expected in the strong held 
regime. 


III. ANALYSIS OF THE DISPERSION RELATION 


In order to gain a better understanding of the dispersion relation © we introduce 
the dimensionless variables and parameters 0 = co/uj cd , K = kc/u pd , = uj pd /uj cdl 

a = (4a/4:57T)(cB 0 /Es) 2 , and S = E 0 /cB 0 . The dispersion relation © then reads 
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where the relativistic gamma factor can be found from the relation 
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We may now discern between the following cases: 

(a) ~ 1: If we let a —> 0, i.e. neglect the QED correction, we only have mode 

propagation for — 1 < 0. When a ^ 0, this is no longer a necessary condition for mode 
propagation, and we can therefore have D — 1 > 0. A special case here is the strongly 
magnetized dusty plasma, cB 0 3> E 0 . Equation © then yields 
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for which the QED effects act only as a small correction. 

(b) 1: From Eqs. (fTUl) and HU we obtain the dispersion relation 
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In the limit a —> 0, we have a single mode given by 
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When a ^ 0, we have further completely distinct modes as given by the dispersion relation 
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As we will demonstrate, Eq. (HU has new branches which is not present in the classical 
theory. Moreover, it will be further pronounced for low frequency waves, and we therefore 
expect the QED effects to play an important role in the propagation of such modes. In the 
limit of large 5, i.e. weak external magnetic held I?o, we can simplify the dispersion relation 
to 
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in order to further elucidate the character of the new QED mode. In the limit a 
obtain the whistler-like dispersion relation D = — K 2 (see Fig. 1). 


0, we 
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IV. APPLICATIONS AND CONCLUSIONS 


In astrophysical environments, such as supernovae, it is believed that charged dust grains 


can affect the plasma dynamics in crucial ways 


It has also been argued that 


supernovae can be a significant source of dust particles |23j], although this claim is under 


debate 


24j . At the same time the plasma may be highly magnetized, thus altering the 


behavior of the dusty plasma significantly 


In fact, the fields can reach amplitudes such 


that the vacuum will be excited, making it nonlinear due to QED effects. In this perspective, 
it is thus natural that the mode structure presented here is more complex than that of linear 
theory. In Fig. 1, we have plotted the dispersion curves, as given by Eq. (USD, for different 
values of the QED parameter A. We note that for K larger than a critical value (which 
depends on A), there are three branches of the dispersion relation dm . in contrast to the 
single branch present in the classical limit. It should be stressed that the values of A used 
in Fig. 1 by no means require unrealistic astrophysical conditions, as the following example 
will show: In interstellar clouds, the dust density is expected to be nod 10 1 m 3 j2a], 
with a typical dust mass m 0 d 10 nTp ro ton fN .J 10 21 kg jjj. For the circularly polarized wave, 
we take the amplitude to satisfy E 0 /Es ~ 2 x 10~ 5 , a value appropriate at a distance of 
1AU from a supernova. In this case, a value Bq ~ 10~ 8 G (well within the limits for such 
environments) results in A ~ 1. For higher dust mass densities we note that even lower wave 
intensities will give values of A close to one. Such dense dusty plasmas are likely to be found 
closer to the core of a supernova remnant. It should also be noted that we have neglected 
radiative-resonant interactions on our new electromagnetic wave mode. However, in many 
cases such interactions can dominate jz/l. fl ], thus altering the behavior of the mode. This 
is a phenomenon that has to be considered in future research. In addition, an important 
problem of the physics of supernova explosions is the relation between the layer of dust 
behind the supernova shock and the process of dust-wave condensation behind the shock 


wave front |22|. Our new nonlinear wave mode may thus play a role from the viewpoint of 
supernova shock propagation. 

We have presented an analysis regarding the propagation of intense circularly polarized 
waves in a magnetized dusty plasma. It is found that corrections due to QED modify the 
mode structure, giving rise to novel low frequency waves with no classical counterpart. Our 
main result is the dispersion relation Q, in which the effects of charged dust and QED are 
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taken into account. The result should be of interest in astrophysical environments, such as 
supernova remnants and interstellar clouds. 
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FIG. 1: The normalized frequency II plotted as a function of the normalized wavenumber K, as 
given by the relation m- Variation of the parameter A = a5 2 Q 2 gives the different solution 
curves. The curves in the upper and lower half planes represents the roots of Eq. (ESP - From the 
line 17 = 0 outwards, the curves for the values A = 0.1,1, 2, 3, and 4 have been plotted. We note 
that the modes for which II > 0 (right-hand polarization) have no classical counterpart. When 
A —> 0 we obtain 17 = — K 2 , which represents a classical limit curve in the lower half plane. 
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